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' O ■ Abstract 

q ■ We present a detailed comparison of the structural predictions of the 

' effective many-body Brenner potential with those of ab-initio studies for 

known reconstructions of diamond (100) and (111) surfaces. These results 
suggest high reliability of the Brenner potential for dealing with carbon- 
based structures where different types of bonding are present at the same 
^v^j , time. 
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1 Introduction 



Carbon-based structures are of a great current interest Q-^]- The challenge of 
these systems from the fundamental point of view is related to the interplay 
between different types of atom bonding, leading to a uniquely large variety of 
structural phases formed by a single element, such as diamond and graphite, 
single- and multi-shell fullcrcnes and nanotubes JjJ and other structures || 
with many peculiar properties. To study the elastic properties and fracture of 
these structures and their mixtures, the transformation paths between them 
[pl§], etc. it is important to develop predictive schemes based on simplified 
empirical potentials, which allow large-scale simulations of complex structures 
with mixed atomic bonding, which are often beyond the possibilities of ab-initio 
calculations. Effective many-body empirical potentials have proven to be useful 
and predictive for a number of materials The potentials developed by 

Tersoff § for group IV elements are very accurate for Si and Ge, also as far as 
interface properties are concerned, but less reliable for C. 

Carbon is particularly difficult for an empirical scheme due to the large 
variety of different types of C-C bonding with very different energetics and 
bond lengths dec- For example, for a single C-C bond in diamond c?cc = 1-54 
A, for a conjugated bond in graphite dec — 1.42 A and for a double bond 
in H2C=CH2 dec = 1-34 A. The Tersoff potential, which has been fit to the 
bulk properties of both diamond and graphite, does not, however, distinguish 
the chemical character of the bond. At diamond surfaces, different types of 
bonding are present at the same time, leading to poor results of the Tersoff 
potential for the surface reconstructions as we show in detail in this work. 

Brenner |Q has re-parametrised the Tersoff potential and added nonlocal 
terms to properly account for the bond modifications induced by a change of 
bonding of neighbouring atoms. As in the Tersoff scheme, the potential energy 
of the system is written as a sum of effective pair terms for each bond, the 
energetics of which depends on the local environment (bond order of Tersoff) 
and, in addition, on the chemical character of the bond (single, double, triple 
or conjugated) derived by evaluating the number of neighbours for the atoms 
forming the bond and all their nearest neighbours. 

Diamond surfaces are an example of a rather simple system, where the inter- 
play between different types of carbon bonding becomes important. Numerous 
calculations exploiting various ab-initio approaches and extensive experimental 
data are available, making the diamond surfaces an important check point to 
verify the accuracy and predictive power of empirical schemes. The Tersoff po- 
tential for C yields the unreconstructed (111) (1 X 1) surface as the most stable 
against the experimental evidence of the (2 x 1) Pandey chain reconstruction Q 
analogous to that of Si(lll). For the (001) face it strongly favours an asymmet- 
ric re-arrangement of carbon atoms beneath the raw of unbuckled dimers || . In 
the present work we compare in detail the predictions of the Brenner and Tersoff 
potentials with results of ab-initio calculations for known reconstructions of the 
diamond(lOO) and (111) surfaces. Our results reveal high quantitative accuracy 
of the Brenner potential. Since the parameters are fit to the bulk properties of 
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Figure 1: Unit cell of (left) the symmetric (2x1) and (right) the asymmetric 
(2xl)a reconstructions of diamond(lOO) after MC energy minimisation with the 
Brenner potential. The top and bottom panels show the top (100) and the side 
(Oil) views. The numbers on the atoms denote the atomic layer they belong 
to. Incquivalent bonds between two atomic layers are labelled with letters a, b 
and c. 



diamond and graphite and to properties of various hydrocarbon molecules, the 
high accuracy at the surface suggests a high predictive power of the potential 
at short distances. With further modifications to include also long-range inter- 
actions (beyond 2 A, the cut-off of the potential) jl0],[ll|], which are now under 
development Jl^j , the Brenner potential promises to become a powerful tool to 
investigate carbon-based structures on a large scale. 



2 (100) surface 

Our previous study Q of the diamond(lOO) surface with the Tersoff potential 
has suggested a new reconstruction with a strongly asymmetric rearrangement 
of atoms in deeper layers. These predictions have now been verified using the 
off-lattice Monte Carlo (MC) technique with the Brenner potential. We have 
confirmed that both the symmetric (2 x 1) and asymmetric (2 x l)a (where 'a' 
stands for 'asymmetric') structures shown in Fig. [I] correspond to local energy 
minima. However, the energies for the two reconstructions given by the Brenner 
potential are found in a reversed order compared to the prediction of the Tersoff 
potential. Relative to the energy of the relaxed (1 x 1) ideal diamond (100) 
surface, the energy gain is found with the Brenner potential to be 5.40 and 4.19 
eV per surface dimer for the (2 x 1) and (2 x l)a structures, respectively. For 
comparison, the Tersoff predictions were 0.26 and 1.55 eV, respectively [||. 

The length of the bonds between atoms in the top four layers are given for 
the two structures in Table l]. For comparison, we also give the results obtained 
with the Tersoff potential |] and those found in ab-initio calculations fl3|| . Note 



3 



Table 1: Bond lengths for the dimerised symmetric (2 x 1) and asymmetric 
(2 x l)a reconstructions (given in A). The bond between an atom in the Nth 
layer and one in the Mth layer is labelled by NM. Whenever more than one 
nonequivalent bond is present, they are labelled by a, 6 and c in Fig. [l] and are 
displayed in the table in this order. For completeness, we show the bond length 
as given by ab-initio calculations [jl3|, the Brenner [this work] and Tersoff Q 
empirical potcnt^s as well as the corresponding bulk diamond bond lcngjths. 



bond 


ab-initio [131 


Brenner [this work] 


Tersoff | E 


(2 x I) 1 " 1 


(2 x 1) 


(2 x l)a 


(2 X 1) 


J2 x l)a 


11 
12 

23 
34 


1.37 
1.50 

1.55, 1.57 

1.56, 1.50 


1.3807 

1.5096 
1.5059, 1.5905 
1.5670, 1.5140 


1.4370 
1.4937, 1.4318 
1.5339, 1.5287, 1.5321 
1.5856, 1.4739 


1.542 

1.515 
1.524, 1.579 
1.570, 1.521 


1.487 
1.496, 1.565 
1.541, 1.489, 1.506 
1.582, 1.477 


bulk 


1.53 


1.5407 


1.5445 



that, from the chemical point of view, each surface atom at the bulk-terminated 
surface has two un-paired electrons (two dangling bonds). Therefore, the dimcr 
bond (bond 11 in Table |]) in the symmetric (2 x 1) surface structure has the 
character of a double-bond. The Brenner potential correctly reproduces the 
length of the bond 11, which is much shorter than both the single C-C bond in 
diamond and the conjugated bond in graphite, but rather close to the length of 
a double bond. It quantitatively agrees with the dimer bond length of 1.37 A for 
the (2 x 1) diamond(lOO) surface found in ab-initio calculations p3|-p^[. Con- 
versely, the Tersoff potential, which does not include the nonlocal terms, predicts 
very different length of the dimer bond; it also gives much smaller reconstruc- 
tion energy for both structures since the chemical character of stronger double 
and conjugated bonds is not accounted for. With the Brenner potential, in the 
asymmetric (2 x l)a structure, the bond 11 is elongated up to 1.437 A (close 
to the graphite value) since it becomes a member of a conjugated (7r-bonded) 
system. Note that the bond 12&, which also connects three-fold coordinated 
atoms, has a very similar length. 

Examination of the first two columns in Table |l| reveals the high accuracy of 
the predictions of the Brenner potential for the bond lengths of the symmetric 
(2 x 1) diamond(lOO) reconstruction. Except a slightly higher difference of the 
length of the bonds 23a and 236, our results agree with the ab-initio results jl3f 
within 0.01 A accuracy. Note that the ab-initio approach of Ref. (k| under- 
estimates the bulk bond length by 0.01 A with respect to the experimentally- 
determined value. In contrast to the surprisingly good quantitative agreement 
of the structural data, the reconstruction energy is different. The value of 5.40 
eV/dimer as given by the Brenner potential is found between the ab-initio val- 
ues (3.02 3.36 jyj and 3.52 |lj| eV/dimer) and the reconstruction energy 
given by the semi-empirical SLAB-MINDO scheme (7.86 eV/dimer |l6||). 
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Figure 2: Side (IlO) views of (left) the relaxed (1 x 1) and (right) the Pandey 
(2 x 1) reconstruction of the (111) surface along with the bond lengths shifts 
relative to the bulk value as given by the Brenner potential. 



3 (HI) surface 



We do not consider the Tersoff potential here since it gives the unreconstructed 
(1 x 1) surface as the minimum energy structure. Fig. || shows the relaxed (1 x 1) 
and the (2x1) Pandey chain reconstructions of the diamond(lll) surface as 
given by the Brenner potential. In brief, the most important changes of the bond 
lengths as compared to the bulk value are as follows. For the relaxed (lxl) 
structure: (i) the contraction of the bond within the first bilayer by —3.5% 
agrees with ab-initio values of -3.1% @, -4.0% @ and -4.2% @; (ii) the 
elongation of the bond between the first and second bilayer by +2.0% agrees 
with +2.1% of Ref. |l7j but seems underestimated compared to +8.7% and 
+9% of Rcfs. JlJ, 19 1 . In the (2 x 1) Pandey reconstruction: (iii) the 7r-bonded 
upper chain bond length of 1.437 A (—6.7%) well compares to ab-initio values 
of 1.47 A (-4.4%) @, 1.44 A |c| and 1.43 A (-6.5%) [pi; (iv) the lower chain 
elongation by +1.4% is close to +0.7% and +0.9% |9|; (v) the stretch of 
the bonds between the first and second bilayers by +3.9% and +4.0% seems to 
be underestimated with respect to the ab-initio values of +8.1% |L7|], +8% 




+4.5% and +6.6% |19 . Further comparison with the results of Refs. ||1 
shows that all other bond shifts agree with the ab-initio calculations within 
~1%. Therefore, except a tendency to underestimate the elongation of the bonds 
between the top and second bilayer, our structural results for diamond(lll) 
agree remarkably well with ab-initio predictions ]l7j-p0|]. 

Relative to the energy of the bulk-terminated diamond(lll), we find energy 
gains per lxl unit cell of 0.244 eV for the (1 x 1) structure (cf 0.37 eV [17| 
and 0.57 eV 111) and 1.102 eV for the Pandey reconstruction (cf 0.47 eV jl7| 
and 1.40 eV $L§). 

We note that there is a long-standing debate on the structural and electronic 
properties of the diamond(lll) surface. An important issue is whether this 
surface is metallic or semiconducting. In most calculations jl7],[nj the band of 
surface states is metallic whereas experimentally the highest occupied state is at 
least 0.5 eV below the Fermi level uM. Dimerisation along the 7r-bonded chain 
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could open the surface gap but only one total-energy calculation obtains slightly 
dimerised chains yielding a 0.3 eV gap p3] in the surface band. Experimentally, 
recent X-ray data |22| does not show any dimerisation but favour the (2 x 1) 
reconstruction accompanied by a strong tilt of the 7r-bonded chains, similar 
to the (2 x 1) reconstruction of Si(lll) and Ge(lll). The tilt is however not 
confirmed by theoretical studies Jl7|, p^,|20[|. 

Neither dimerisation nor buckling of the 7r chain is found in our results for 
the Pandey reconstruction in agreement to most ab-initio results. However, our 
recent MC study Q of the structure of diamond(lll) based on the Brenner 
potential has shown that, in addition to the stable (2x1) Pandey chain re- 
construction, there exist additional meta-stable states, specific for carbon, with 
all surface atoms in three-fold graphite-like bonding. Since the energy of these 
metastable states is very close to the one of the Pandey (2 x 1), these struc- 
tures can coexist with the Pandey structure at a real surface. Moreover, due to 
symmetry breaking induced by a strong dimerisation of the lower (4-fold coordi- 
nated) atomic chain in the first bilayer, the meta-stable reconstructions is likely 
to exhibit semiconducting behaviour. Although the new structures and their 
surface electronic properties ought to be checked in ab-initio studies, the high 
accuracy of the Brenner potential demonstrated in this work strongly supports 
this prediction. 



4 Conclusion 

We have performed off-lattice Monte Carlo study of the (100) and (111) diamond 
surfaces with the empirical many-body Brenner potential |7j] and compared the 
results in detail with those obtained with the Tersoff potential |9| and with 
ab-initio approaches |l3|-|l5|, P^|~|lS|] ■ We find that the Brenner potential is ex- 
tremely accurate in describing the structural properties at surfaces, supporting 
the recent predictions [^3| of new meta-stable reconstructions of diamond(lll). 
On the other hand, the Tersoff potential (H), which does not distinguish the 
chemical character of the bond, turns out to give a poor description of surface 
properties. 

The Brenner potential, however, cannot describe weaker long-range interac- 
tions, such as the interplanar interactions in graphite due to the cut-off at 2 A. 
This is the most serious limitation to be overcome. Recently, further modifica- 
tions of the Brenner potential to include also long-range interactions beyond 2 
A are being proposed 10-1^]. Given the high accuracy of the short-range part, 
the modified Brenner potential IO-I4] promises to become a powerful tool to 
investigate carbon-based structures on a large scale. 
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